Nuclear Symmetry Energy from QCD sum rules 
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We calculate the nucleon self-energies in an isospin asymmetric nuclear matter using QCD sum 
rules. Taking the difference of these for the neutron and proton enables us to express the potential 
part of the nuclear symmetry energy in terms of local operators. We find that the scalar (vector) 
self energy part gives negative (positive) contribution to the nuclear symmetry energy which are 
consistent with the results from relativistic mean field theories. Moreover, we find that an important 
contribution to the negative contribution of the scalar self energy comes from the twist-4 matrix 
elements, whose leading density dependence can be extracted from deep inelastic scattering experi- 
ments. This suggests that twist-4 contribution partly mimics the exchange of the <5 meson and that 
it constitutes an essential part in the origin of the nuclear symmetry energy from QCD. Our result 
also extends an early success of QCD sum rule method in understanding the symmetric nuclear 
matter in terms of QCD variables to the asymmetric nuclear matter case. 
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I. INTRODUCTION 

With the recent constructions and plans for the next 
generation low energy rare isotope accelerators world- 
wide, there is a renewed interest in the nuclear symmetry 
energy, as the topic could be studied further with finer 
details in these experiments Understanding the de- 
tails of the nuclear symmetry energy will provide valuable 
insights into exotic nuclear matter ranging from rare iso- 
topes to neutron rich nuclear matter such as the neutron 
star 0, Q ■ One of the main puzzle to be solved is the 
behavior of the nuclear symmetry energy at high den- 
sity gi. 

From a phcnomenological point of view, insights into 
nuclear symmetry energy can be obtained from under- 
standing the nuclear binding expressed through the semi 
empirical mass formula [f| . There, the symmetry energy 
can be understood as originating from the energy differ- 
ence between the proton and the neutron in an isospin 
asymmetric nuclear matter. Hence, within the approxi- 
mation, the problem of calculating the nuclear symmetry 
energy can be reduced to obtaining the optical potential 
and thereby calculating the energy of the nucleon quasi- 
particle near the Fermi surfaces in an asymmetric nuclear 
matter. 

In Dirac phenomenology of nucleon-nucleus scatter- 
ing @, the optical potential of the nucleon is com- 
posed of a vector and scalar part, U ~ S + Vj°. It is 
well known that in order to fit the spin observables, one 
needs a strong scalar attraction (Re S < 0) and a strong 
vector repulsion (Re V > 0) both of several hundred 
MeV, but such that the combined sum gives only a small 
attraction of few tens of MeV, consistent with the tradi- 
tional low energy nuclear physics. The strong scalar and 
vector potential comes about naturally also in relativis- 
tic mean field theories (RMFT), where meson exchange 
interaction between nucleons on the Fermi sea produce 
strong scalar and vector potentials for the nucleons. 

But it was only after the works in QCD sum rules that 



the strong optical potentials were found to have a basis 
from QCD. The first pioneering work of applying QCD 
sum rule method to nucleons in medium were performed 
by Drukarev and Levin [t| [l(j. Later, the relation be- 
came clearer through the work by Cohen, Furnstahl and 
Griegel (Tlj . who used the energy dispersion relation to 
show that the showed the strong scalar- vector self energy 
appearing in the quasi-nucleon pole in the symmetric nu- 
clear matter jn] - [l4j | can be traced back to the scalar- 
vector quark condensate in the nuclear medium and that 
the contribution from higher dimensional condensates do 
not change the general behavior. 

Motivated by these results, and to express and eluci- 
date the origin of nuclear symmetry energy directly from 
QCD, we have applied QCD sum rule to calculate the 
neutron and proton energy in the asymmetric nuclear 
matter. Identifying the difference with appropriate fac- 
tors to the nuclear symmetry energy, we show that this 
energy can be expressed in terms of quark and gluon de- 
gree of freedom. Attempts to calculate the nucleon mass 
in an asymmetric matter using QCD sum rules were re- 
ported before [TH - flTl ]. But here, we will follow the for- 
malism adopted in Ref. [lljj. We have performed the 
operator product expansion (OPE) up to dimension 6 
operators and have identified all the independent twist-4 
operators, where for most of the operators, the leading 
density dependence can be extracted from the deep in- 
elastic scattering data on higher twist effects. From the 
QCD sum rule analysis, we find that the scalar (vector) 
self energy part gives negative (positive) contribution to 
the nuclear symmetry energy which are consistent with 
the results from relativistic mean field theories. More- 
over, we find that an important contribution to the neg- 
ative contribution of the scalar self energy comes from 
the twist-4 matrix elements, whose higher density be- 
havior will determine the still controversial property of 
the symmetry energy at these densities. Our result also 
extends an early success of QCD sum rule method in 
understanding the symmetric nuclear matter in terms of 
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QCD operators to the asymmetric nuclear matter case. 

The paper is organized as follows. In Sec. II, we start 
with brief review and simple idea for the nuclear symme- 
try energy. QCD sum rule for nucleons in the asymmetric 
nuclear matter is presented in Sec. III. Finally, sum rule 
analysis and conclusion are given in Sec. IV. 



by from the second term Eq. ([2]). That is, (E n — E p ) = 
\IAA, hence, 

a A =±AA = ± } (E n (AJ)-E p (AJ)). (4) 

For infinite nuclear matter case, one can get Epj as: 



II. A SIMPLIFIED DESCRIPTION FOR THE 
NUCLEAR SYMMETRY ENERGY 
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N 



J" d^ kp 



d 3 k n d 3 k p E N (p n , p p ). 



(5) 



We start from a finite nuclei with A nucleons. In the 
Bethe-Weizsaker formula for the nuclear binding energy 
given as 



Zm r 



E B /c 2 , 



mtot = Nm r _ 
E B = a v A - a s A^ - a c {Z{Z 
-a A I 2 A + 6(A,Z), 



1))A- 



(1) 



where / = (N — Z)/A, the fourth term accounts for the 
total shifted energy due to the neutron number excess. 
Taking the infinite nuclear matter limit of this formula, 
one notes that a A reduces to the nuclear symmetry en- 
ergy [I]. 

To derive the formula for a A that can be generalized 
to the infinite nuclear matter, we start from a simple 
formula for the total energy: 
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A{\ + I)E n - 

A(E n + E p ) 



1 



A(l-I)E P 
-Al(E n -E p ), 



(2) 



where E n (E p ) is the average neutron (proton) quasi par- 
ticle energy in the asymmetric nuclear matter. Now, the 
core of the model is what approximation goes into calcu- 
lating the average energy. 

The symmetry energy in an asymmetric nuclear matter 
is defined as 



And the nuclear symmetry energy can be obtained by col- 
lecting coefficients of I 2 in J2]). This generalized E sym (p) 
in Eq. ((3]) can be decomposed into kinetic like part and 
potential like part by mean field type quasi-particle ap- 
proximation, in which one can treat both parts sepa- 
rately. The kinetic part of E sym can be obtained ac- 
cording to Ref. [l8j : 



k 2 F + E^ V ^ I=Q ^ 



(0) 



where kp is the Fermi momentum and E q V ^ I= ^ is the 
potential part of the quasi-nucleon self energy in sym- 
metric nuclear matter at saturation density. 



A. Linear density approximation 

In QCD sum rule calculations, we will be using the lin- 
ear density approximation, because in medium conden- 
sates in QCD sum rule can be most reliably estimated to 
leading order in density. This means that for either the 
proton or the neutron, the mass will be given as follows: 

Eyipn, P P ) = m + ap p + bp n 

= m o + 7}P(a + b) + ^pl{b - a), 
E v(Pn,P P ) =m + ]-p{a + b) - \pl{b - a), (7) 



E tot (p,I)^E ( P )A + E syni (p)I 2 A + O(I i ), (3) 

where p is the nuclear medium density and / = (N — 
Z)/A — > (p n — P P )/{pn + Pp) and the neutron and proton 
densities are p n = + 1), Pp = \p(\ — -0 respectively. 
Therefore, in Eq. ([2]), the symmetry energy will have con- 
tributions from the term proportional to / in (E n — E p ) 
and the term proportional to I 2 in (E n + E p ). 

For a non interacting Fermi gas of nucleons with mass 
ton, calculating the average nucleon energy will give E = 
^Ep, where Ep is the nucleon Fermi energy. Following 
the procedure described above and extracting the term 
proportional to I 2 gives a nuclear symmetry energy of 
\E F .^ 

Going back to finite nuclei, assuming a 'Fermi well' 
with constant energy difference A between adjacent nu- 
cleon energy level, the symmetry energy can be obtained 



where too is the vacuum mass and a, b are the constants 
to be determined later. We then have, 



E 



■N 



1 



J" d^ d^ kp 

i 



i 



m o + ^Pp + 2 bf)7 



toq + -p(a + b)+ -pl(b - a) 



Combining Eq. ([5]) with Eq. ©, we obtain the sym- 
metry energy. That is, (E^ -E^) = \{E^(p n ,p p ) - 
E v(Pn,P P )), hence, 



ET = -(EZ( Pn ,P P ) ~ E p v (p n ,p p )), (9) 
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is similar to the relation given in Eq. (|4]). Therefore, 
to this order, the symmetry energy comes only from 
the energy difference in the proton and neutron at the 
Fermi surface in an asymmetric nuclear matter as given 
in Eq. ([9]). However, for operators with higher dimension 
and higher density dependence, the factors appearing in 
Eq. ([SJ should be modified and symmetry energy will 
have contributions from both the sum and the difference 
of the nucleon energies. 

The quantity of interest, namely, (E v (p n ,p p ) — 
Ey(p ni p p )) can be obtained by looking at the pole of 
the nucleon propagator in nuclear medium: 

G(q) = -i [ d 4 xe^(vI/ |T^(x)^(0)]|vI/o), (10) 



where I'J'o) is the nuclear medium ground state, and ip(x) 
is a nucleon field. Relativistic mean field type of contri- 
bution will then appear in the self energies. The nucleon 
propagator can be decomposed as 



G(q) = G s (q 2 ,q ■ u) + G q (q 2 , q ■ u)4 + G u (q 2 ,q ■ u)i, 



III. QCD SUM RULE FOR NUCLEONS IN THE 
ASYMMETRIC NUCLEAR MEDIUM 



A. Operator Product Expansion and Borel sum 
rule 

To express the self energies in terms of QCD variables, 
we start with analyzing the correlation function via the 
operator product expansion (OPE). The Correlator is de- 
fined as 



n(g) 



d A xe igx < 



* o \T[r)(x)f}(0)]\*o), (18) 



(11) 



where rj{x) is an interpolating field of nucleon and l^o) is 
the ground state of the asymmetric nuclear medium la- 
beled with the rest frame medium density p, the medium 
four-velocity and the asymmetric factor I. This 
ground state is assumed to be invariant under parity and 
time reversal. We will be using the Ioffc nucleon interpo- 
lating field given as [ill EH > 



where is the four-velocity of the nuclear medium 
ground state [l2j . 

The nucleon self energy can be decomposed similarly 
as FMi: 



H{<l) = ^s{q\q-u) + t»{q) li 



(12) 



r)(x) = e a bc[ul (x)C7 M u b (a;)]757 M (i c (a;). 



(19) 



As in the case of the nucleon propagator, using Lorcntz 
covariance, parity and time reversal, one can decompose 
the correlator into three invariants [131 : 



where 



E{f(?) = £ u (g 2 , q ■ u)u" + £ g (g 2 , q ■ u)q^ 



(13) 



IL(q) = n s (g 2 , q ■ u) + U q {q 2 ,q ■ u)4+n u (q 2 , q ■ u) 



(20) 



In the mean field approximation E s and £„ are real and 
momentum independent, and T, q is negligible. Hence, 



1-S a 



1 - So 



(14) 



The phenomcnological representation of the nucleon 
propagator can then be written as 



G(q) = 



1 



4-M n -Y.(q) 



A 2 i+M 



(go - E q )(q - Eq) 



7=TT> (15) 



The three invariants are functions of q 2 and q ■ u, while 
vacuum invariants depends only on q 2 . For convenience, 
we set the rest frame for the nuclear medium, which 
means — > (1,0); also we fix \q\. Hi(q 2 ,q ■ u) then be- 
comes a function of qo only, which means n^g 2 , q ■ u) — > 
IIi(go, |g| fixed) (t = {s,q,u}). 

As mentioned before, we will follow the formalism 
adopted in Ref. [ll[ and write the energy dispersion rela- 
tion for these invariant functions at fixed three momen- 
tum \q\: 



1 f°° 

A is unity in this discussion, but if one did not neglect ^U\Q0: = J duj 



H q , effects of T, q should be accounted in A 2 as the factor 
(1 — E ) _1 . E q and E q are the positive and negative 
energy pole: 



w - g 



polynomials, 



Eg =Z V + X /<f+Mtf, 



Eq — 



f+M* N 2 . 



(16) 
(17) 



An,(w, \q\) = lim [Hi(cj + ie, \q\) - U,(tu - ie, 

£->0 + 

= 2Im[n i (o; ) |g1)]. 



(21) 
(22) 



With fixed |gj, G(q) only depends on g . One can extract 
self energy near ~ E q with analytic properties of the 
nucleon propagator. 



The lowest energy contribution to the discontinuity 
will be saturated by a quasi-nucleon and quasi-hole con- 
tribution in the positive and negative energy domain re- 
spectively. Their contribution to the spectral density will 
be given as in Eq. (|15[) . which will have the following con- 
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tribution to the invariant functions: 

n s (g ,|g1) 



x* 2 



H q (q ,\q\) = -X* N 2 

n„(g 0! k1) = +A^ 



M*N 


(<?o 


-E q )(q 


-E q ) 




1 




(<7o 


-E q )(q 


~E q ) 








(<?0 


-E q )(q 


~E q ) 



+ 



(23) 
(24) 
(25) 



where A J 2 is the residue at the quasi-nucleon pole, which 
accounts for the coupling of the interpolating field to the 
quasi-nucleon excitation state, and the omitted parts are 
the contributions from the higher excitation states, which 
will be accounted for through the continuum contribution 
after the Borel transformation. 

The even and odd parts of the invariant functions are 
respectively related to the following parts of the discon- 
tinuity: 



W(Qo,\Q\) 



2ni 



cllu 5 ^ h polynomials, 

u 2 - Qo 

duj ^ — ' 5 + polynomials, 



^ 2 ~Qo 



(26) 



where we have defined invariants as follow for conve- 
nience 



Il^go, |<?1) = nf (q 2 , |g1) + q I^ \q 2 , \q\). (27) 



The OPE of the three invariants of both the even and 
odd parts can be expressed as 



Ili(q 2 ,q 2 )=Y / Ci l (q 2 ,q 2 )(O n } p j, 



(28) 



where {O n ) p j is the ground state expectation value of 
the physical operator in the asymmetric nuclear medium; 
(*o|O n |* )pj- We will be adopting the OPE at q 2 -)• 
— oo at finite |g| — > fixed; this is equivalent to the limit 
of g 2 —> — oo at finite |g| — > fixed. The Wilson coeffi- 
cients CI (q 2 , go) can thus be calculated in QCD at short 
time [n(. 

The OPE of the correlator for the proton interpolating 
field up to dimension 5 operators are given as follows: 



nf(g 2 ,k1) 
n s °(<? 2 ,k1) 



nf(g 2 ,k1) 



-^g 2 \n(-q 2 ){dd) pJ + ^%(d{iD iD }d) p ,i, 
47r z 6ir z q z 



— \n(-q 2 )(diD Q d) p j, 
(g 2 ) 2 ln(-g 2 ) 



64tt 4 
( 1 



+ (g^Hz 2 ) - if) (d{^D }d) pJ + (^H-q 2 ) ^f) (u{^D }u) pJ 

1 W_ga) (^±G 2 ) - ( \n(-q 2 ) 

Itt v y ; \ 7T I p,I 144tt 2 V v y ' 



4g 2 \ /a 



32 



n ?(% 2 , kl) = H-q 2 ) [(^u) p j + (dU) 



p [(M . G) 2 + (u . G) 2 ] 

q- I \ 7T / p ,i 



2 q 2 

J - 3n 2 (g 2)2 (uboiD iD }u) p ,i 



2 q 2 2 1 11 

-^J^( d {^ iL> o iD o}d) P j - ^2 ^2 («{7o*A)«'-Do + (gsu^a ■ Qu) p>1 , 



(29) 
(30) 



-Lg 2 ln(-g 2 ) [7{u^u) pJ + (dU) pJ ] + -| % {u{ lQ iD Q iD }u) p ,! + ^%<J{ 7 oi£>oiA)}d> Pl / 

127T 2 J tt z q 1 tt z q z 

7T2 ln (-9 2 )(.9suV • Qu) p j + — -Ur In (- g 2 ) (,9s tf T a • 

f>ir^ 127T^ 



n?(g 2 , 151) = - g^ln(- g 2 )(d{ 7 oi£>o}d) p ,/ - ^ln(-g 2 )(w{7o» J Do}^)p,/ 



+ 1 l n( V)(-[(^G) 2 + (u.G) 2 ] 

307T Z \ 7T 



I 



(31) 



(32) 



(33) 



(34) 



The quark part and their flavor structure of the above OPE can be obtained by suitable substitutions of the cor- 
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responding OPE for the S given in Ref. [20( ; by changing 
q — > u, s — > d, and neglecting terms proportional to m s . 
Moreover, when both u and d quarks arc identified to the 
generic light flavor q, our OPE also reduces to that given 
in Ref. Q- 

Next task is identifying the nucleon self energies in the 
asymmetric nuclear medium. So we have to concentrate 
on the quasi-nucleon pole and not on the quasi-hole nor 
on the continuum excitations. To this end, we apply the 
Borel transformation with ap pro priate weighting func- 
tion to the dispersion relation |l2|: 

B[n t ( qo , \q\)} = — du Wi^AIUiu, \q\), (35) 
2™ J_ Uo 

(36) 

where E q is the quasi-hole pole which will be assigned 



W{u) = {u-E q )e- J> ' M *, 



as phcnomcnological input. The weighting function will 
de-emphasize the contribution from the quasi-hole, and 
the Borel transformation suppress the continuum contri- 
bution. Using Eq. (|26p . the OPE side of the sum rule 
can be obtained by taking the Borel transformation of 
11* (go, \q\) = nf (q 2 , \q\) - E q X\° (q 2 , \<f\). Here, we define 
differential operator B for the Borel transformation of the 
OPE side as 



B[f{qlM)] 



lim 

— <2q ,n— >oo 
-q 2 Jn=M 2 

f(M 2 ,\q\), 



(-9o 2 )" +1 



d_ 

dq'o 



f(ql\Q\) 



(37) 



where M is the Borel mass [21[ . 

The Borel transformed invariants which contain con- 
tinuum corrections arc as follows: 



B[U s (ql\q])} = 



y 2 M; e -^-^ M2 



4tt 2 
E„ 



(M 2 ) 2 Ei(dd) pJ 



4 
3^ 



2 f(d{iD iD }d) p ,iL-° 



;M 2 E (diD d) p jL- 



(38) 



23[II g (g 2 ,|g1)] = 



{M 2 fE 2 L~ 



2,2-k 1 



1 

9^ 



M 2 E a - —q 2 ) (d{ 10 iDo}d) pJ L-i 



^ M ' 2E ° ~ g-2? ) {u{wDo}v) pJ L-i 



1 

327T^ \ 7T I p.I 

1 



M 2 (^G 2 ) t E L-* + (M 2 Eo - 4?) ( ^[(« ■ G) 2 + (u ■ G) 2 } ) L' 



p,i 



E n 



6tt 2 



M 2 E Q L-* [(u^u) p j + (dU) p j] -JLfl-JLf) 



^2 f 1 - |^) (di-roiDoiDoWpjL-i 

2 1 

— (u{'~foiDoiDo}u}pjL~» + -^-^(g s v)(r ■ Qu) p jL~~> 



(39) 



B[IV u {ql 



= X* 2 ^e-( E ' 2 -^ M2 



1 3 

[7(v?u)p,i + {d f d) p j] E-Jrs + —f {u{ l0 iD iD Q }u) pjL-i 



12tt 2 (M 2 ) 



^(f{d{ l0 iD iD Q }d)pjL-r> - —M 2 (g s u^a ■ gu) PtI E L- 



—rM 2 {g s d^a-gd)pjE L- 



E n 



16 



- ^M' 2 (d{ lQ iD Q }d)pjE L-t - ^M 2 (u{ lQ iD Q }u)pjE Q L- 
1 



(40) 
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Here, we include the contributions from the anomalous- 
dimensions as 



-2r„ 



-r 



HM/Kqcd) 



\ii(h/Aqcd) 



-2r, ; +r „ 



(41) 



where T„ (To n ) is the anomalous dimension of the inter- 
polating held i] (O n ), p is the normalization point of the 
OPE, and Kqcd is the QCD scale [HQ. 

Also, the continuum corrections are taken into account 
through the factors 



E Q = l-e s °/ M '', 

Et = 1 - e s *>/ M ' 2 (s*/M 2 + 1) , 



E-, 



(42) 
(43) 
(44) 



where Sq = 



Q—q 2 , and ujq is the energy at the continuum 
threshold. As in many previous studies, we choose ujq = 
1.5 GeV. 



I 

B. Condensates in the asymmetric nuclear medium 



To estimate the matrix elements, we will use the linear density approximation in the asymmetric nuclear matter: 

(6) pJ = (0) vac + {n\6\n)p n + (p\d\p) Pp 

= (O)vac + \ ((n\6\n) + (p|d|p>) p + \ ((n\0\n) - (p|6|p>) I p. (45) 

The quark flavor of condensate becomes important in the asymmetric nuclear medium. Consider an operator O u ^d 
composed of either u or d quark respectively. Making use of the isospin symmetry relation, 

(n\6 Utd \n) = (p\6 djU \p), (46) 

we can convert the neutron expectation value to the proton expectation value, thereby rewriting Eq. (|45p for the two 
quark operators as follows: 

(duM) P j = {d u4 ) v&c + ({p\d Q \p) t (p|6ib)i>. (47) 

Here, ' — ' and '+' are for and 'd' quark flavor respectively, and the isospin operators are defined as 

do = ^(d u + 6 d ), 6 1 = ^(6 u -d d ). (48) 

Hence, we will convert all the expectation values in terms of the proton counterparts and denote them as (p\0\p) — > 
(0) p , throughout this paper. The next task is to find (Oq) p and {0\) p for all operators appearing in our OPE. 



1. {qD^ 1 ■ ■ ■ Dfj, n q) type of condensates 

Let us start by estimating the lowest dimensional operators ([qq]a) P and ([qq]i) p . To find ([qq]i) p , we will use an 
estimate based on using the QCD energy momentum tensor in the baryon octet mass relation to leading order in the 
quark mass[22j]; Eq. (|A3|) in Appendix A. Using Eq. (|A4[) , one finds 



mo, - \ ((**> - wm - \ ^ + l„!^ +m ^ ) ■ w 

We will use the baryon masses as given in the Particle Data Group |23j]: m^o = 1315 MeV, mg- = 1321 MeV, 
m s+ = 1190 MeV, m s - = 1197 MeV. Using m s = 150 MeV and m q = \ (m u + m d ) = 5 MeV, Eq. (gSJ) becomes 
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For {[qq]o)p, we make use of the nucleon <jn = 45 McV term: 

(Mo) P = \ {(p\uu\p) + (p\dd\p)) = £L ~ 4.5. (51) 

For convenience, one can introduce the parameter lZ±(m q ) defined as 

{p\uu\p) ± (p\dd\p) = U±(m q ) ■ (p\uu\p), (52) 

which leads to 

Using the previously selected values with the explicit quark mass dependence, we have 

/ ,((tn 249 MeV \ / f a N 249 MeV \\ 

K±{m q ) = \1 ± ^— - 30Q _ 2m J / [— + 300 MeV _ 2m J J > i 54 ) 

so that 7£±(m 9 = 5 McV)) = 1 ± 0.68. 

Using this parametrization, we can express the u quark or d quark condensates as follows, 

([qq}u.d)pj = ([5<?]tt,d)vac + U =F ^j^j 7 ) ' ([w]o) P (55) 

where [qq] u = uu and [qq]d = dd. For (qq) vac , wc use the Gellmann-Oakcs-Renner relation: 

2m q (qq) vac = -mlf%, (56) 

where = 138 McV and f n = 98 MeV [3. For m q = 5 MeV, we have (gg) V ac = -(263 McV) 3 

Likewise, we will further assume that the ratios between the isospin singlet and triplet operators remain the same 
for all two quark operator expectation values with any number of covariant derivatives inserted: 

{[qD^-.-D^q]^ = ■ {[qD^ ■ ■ ■ D^q] ) p . (57) 



With this assumption, {qD pl ■ ■ ■ D fln q) p j can be written as 



U-(m q) 



{[qD pi ■ ■ ■ D^q] uA ) p j = {[qD M1 ■ ■ ■ D Pn q] uA ) vac + \ l =p n+ ^ m j 1 J ' ([l D m ' ' ' D ^nQ\o)p P- (58) 

The symmetric and traceless part of the above type of expectation values, constitute the moments of the twist-3 
e„(x,/i 2 ) structure function defined as follows [36| : 

([q{D fll ---D^}q} ) p = (-i) n e n ( f x 2 ){p fll ---p fln }, (59) 

e n (n 2 ) = dx x n e n (x, fi 2 ), (60) 



where {fi\ ■ ■ ■ fi n } means symmetric and traceless indices. Then the two quark twist-3 condensates in our sum rule 
can be written as follows: 

{[qiD^q} u4 ) p j = {[qiD Q q] Utd ) p ,i • u'^ = m q ([q^ ' q] u ,d) P .i = 0, (61) 
{[q{iD p ,iD v ,}q] u . d ) p j = ^ {[q{iD iD Q }q] u>d ) p>1 ■ ^< - |S7«/J , (62) 

where the in- medium rest frame u' = (1,0) has been taken and the matrix element is estimated as 

{[q{iD p ,iD y i}q] a ) p = M%e 2 {p 2 ) (u' p u' u - jft^J , (63) 



where one can identify that Mf ! e2(p 2 ) = §([<z{*-Do*-Do}<z]o}p an d {[q{iDoiDo}q] Ui d.) p ,i can be written as 

([q{iD iD }q] u4 ) pJ ~ ( 1 T !4^4 7 ) ' M N^ 2 ) P- (64) 

Since there are no measurements on the twist-3 structure function, we will take the estimate for M^e2(p 2 ) ~ 0.3 GcV 2 
given in dHH. 

When spin indices are contracted, the operator becomes 



([qD 2 q} u4 ) pJ = \([g s qcr ■ Gq} u ,d)pJ = \ { 1 T 35— 4-M • ([ffs^ ■ £g]o) P P, 
2 2 \ U+(m q ) J 

where {[g s qcr ■ Gq]o) P is chosen to be 3 GeV 2 as in Ref. [Til. l26j. 



(65) 



& (?7mi -^M2 ' ' ' -^Mn *J/P e °/ condensates 
The simplest condensate of this type is 

(qj\q) P j = (qi'q) P j u '\ -> (qU) P ,iu'\- (66) 

For this, the ratio (u^ u) p / {ct 1 d) p = 2, and the iso-spin relation for (q^q) p ,i can be written as 

([? f «]i)p = ^([<Z + <7]o>p, (67) 
which leads to the following matrix elements appearing in the sum rule: 

(kUU) P ,i = (l T • ([qUh) P P=(^T \lj P- (68) 

When covariant derivatives are included, one can estimate the two quark twist-2 condensates from the corresponding 
parton distribution function: 



(—i)n-l 

(q{l Pl D P2 ■ ■ ■ D Pn }q)v = -J^ A n^P^ ' * '^J, 
where A* {p 2 ) = (A£(fi?) + Af i (p, 2 ))/2 is the reduced matrix element [37|,[38|]: 

AI(P 2 ) = 2 / 1 d«' , - 1 [ g ( I ,^) + (-l) n 9>^ 2 )], (70) 



where g(x,/x 2 ) and q(x,p, 2 ) are the distribution functions for quarks and antiquarks in the proton respectively, /i 2 is 
the renormalization scale. For the distribution functions, we used the leading order (LO) parametrization given in 

Ref. m. 

Specifically, the spin-2 part can be written as 

{[q~{lniD v }q] Ut d) p ,i ->■ {[q{l p 'iD yl }q] u , d ) pJ = ^{[q{y iD }q] Ujd ) p>1 ■ (u' p u' v - ^9 P ^j , (71) 

where the in-medium rest frame has been taken. The matrix elements for each flavor in {[q{ r yoiDo}q\ u _ c i) p can be 
identified as 

{u{ lp aD v ,}u) p = Imjv^Cm 2 ) • K< - \g„u) , (72) 
{d{ ltl ,iD v ,}d) p = ±M N A d 2 (n 2 ) ■ (u'^u' v - l 9ia \ , (73) 
where A^y 2 ) ~ 0.74 and A^p?) ~ 0.36 at p 2 = 0.25 GcV 2 (LO) HH- 
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(76) 



(77) 



One can introduce a ratio factor for (0\) p as 

{[q{-l^iD vl }q] l ) p = UaM) ' <[9{7m'*A4<z]o>p, (74) 
where H A2 (/j 2 ) = (A% - A$)/(A% + A%) ~ 0.35 so that ([g{7o«A)kkd)pj 

can be written as 

([9{7o«A)}4m)pJ = (l ^TZ A2 (p 2 )l) ■ ([q{joiDo}q\ ) p p 

= (lT^^ 2 )/) • hl N Al(p 2 )p. (75) 

The spin- 3 part can be written as 

{[q{l\'iD^iD vl }q] Utd )pj = 2{[q{^ Q iD a iD }q] Utd ) pJ ■ (u x u'^u' u - ^(u^g^ + u'^gxv + u' v gx^)j , 
where the matrix elements for each flavor in ([g{7o*-Doi-Do}<z]ti.d)p can be identified with 

{u{yx'iD^iD v >}u) p = ^M N A%(p?) ■ (u' x u'^u' v - ^(u x g^ + u'^gxu + u' v 9\S) > 
{d{yx'iD^'iDu'}d) p = ^M N Aj(p 2 ) ■ (u' x u'^u' v - ^(u' x g^ + u'^gxu + u' v 9\S) > ( 78 ) 

where A^p 2 ) ~ 0.22 and Ai(p 2 ) ~ 0.07 at p 2 = 0.25 GeV 2 (LO) [H- Similar to the spin- 2 condensate case, one can 
write {0\)p for spin-3 condensate as 

{[q{ 1 X'iD lll iD v ,}q] 1 )p= K A3 (p 2 ) ■ {[q{lx>iD„>iD„,}q\ } p , (79) 

where TIa 3 (h 2 ) = (A% - Af)/(A% + Af) ~ 0.5f and {[q{j iD iD a }q} u , d ) pJ can be written as 

([q{loiD iD }q\ Utd ) pJ = (l qp Ka 3 (p 2 )l) • ([?{7o*A)iA)}<?]o)p p 

= (f =F K A3 {p: 2 )l) ■ l -M N Al{p 2 )p. (80) 

Operators with contracted spin indices are 

([qpq] u ,d) P ,i = 0, (81) 
([q t D 2 q] u . d ) pJ = -([ff a g + <7 • Gq]u,d) P ,i - ^{ l T T^aJ) • ([g s q^<J ■ Qq] Q ) p p, (82) 

where ([g^gV • Gq]o)p is chosen to be —0.33 GeV 2 flil. |26|. 

3. Gluon condensates 

As for the gluon operators, because they do not carry quark flavors, the expectation values do not depend on /. 
These operators can be written as [HI, [l4[ 



ri G *\ =/^G 2 ) -2(^{E 2 -B 2 )) p, (83) 

7T / p,I \ 7T / vac \ 7T ' p 

^[( U .G) 2 + ( U .G) 2 ]) = -(^(£ 2 + £> 2 )) p, (84) 

7T / p,I \ IT I p 

where _B and B are the color electric and color magnetic fields. For the expectation values of the gluon operators 
we take; ((a s /7r)G 2 ) vac = (0.33 GeV) 4 [13, {{a s /Ti){E 2 - B 2 )) p = 0.325 ± 0.075 GeV and ({a s /n){E 2 + B 2 )) p = 
0.10 ±0.01 GeV pf. 
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C. Dimension 6 four-quark operators 

In many previous QCD sum rule studies, dimension 6 four-quark condensates are assumed to have the factorized 
form as 

(u a a u b ^u d s ) pJ ~ (u a a u b p ) pJ (u^u d s ) p j - (u a a u$) p j{ U Cu*} PtI , (85) 
Ku b p d c 7 d d 5 ) pJ ~ <<^} p , 7 <^) pJ . (86) 

While large N c arguments can be made to justify factorization in the vacuum, no such argument exists in the medium. 
Using the following steps, one can classify the four-quark condensates in terms of the independent operators and of 
different twist. 

1. Twist-4 operators with pure quark flavor 

The first type of four-quark operator appearing in the OPE of the nucleon sum rule, involves quark operators with 
the same flavor, and is of the color anti-triplet diquark times triplet anti diquark form. Using the following Fierz 
transformation, one can identify the independent four-quark operators in terms of products of quark-antiquark pairs. 

iabcia'b'c(UaCj p U b )(ul,^„Cu a ') = f-abc^a'b' c — (u a 'T° U a )(u b 'T k U b ) ■ Tr [^ p T k^vCT^ C] 

= tabcta'b'c ^(Ua'U a )(Ub'Ub) ' (-4^) 

+ (u a >j 5 u a )(u b ,~f 5 u b ) ■ {Ag pv ) 
+ (u a ,"f a u a )(u b ,"f f3 u b ) ■ {4S P i3„ a ) 

~ (w a <7 Q 75"a)(Wfc'7' 3 75W6) ' ( 45 W^a) 

+ {uaio aa u a )(u b >(j f3 Pu b ) ■ -Tr [7 Al a (3( g7 I y(T Q a] 

+ (Wa'7" u a)(u&'7 /3 75Wb) ' (8l£ M 
~ {u a >U a ){u b ,(T aa U b ) ■ {&ig afl gav) 

- (u a >l5U a )(u b ,<T aa u b ) ■ (4e /ivaa )|, (87) 

where T = {I, y a , 170,75, cr Q(3 , 75} and = g pa 9vp + 9^9au ~ 9^9a/3- 

When quarks of the same flavor combine into a diquark, certain combinations are not allowed due to Fermi statis- 
tics. From these conditions, one can extract constraints among four-quark operators that can be used to identify 
independent operators. Among several conditions, the most suitable constraint for our OPE can be obtained from 
the zero identity used in Ref. [35|. With the constraint Eq. (|B2[) in Appendix B, Eq. (|57|) can be simplified as 

£abcta'b>c(uTC^ p U b )(ul,^ u Cu a i) = e a bc^a'b'c \ [(u a ^ a U a ) {u b , -J P U b ) - (u a ^ a -/ 5 U a ) {u b , ^ ^U b )] ■ (8S paua ) 



+ (w a '7 a '«a)(u&'7 /3 75'«&) ' (leie^a) 



The last term in Eq. (|88j) will be dropped as one should take expectation value with respect to parity-even nuclear 
medium ground state. Then, only two types of four-quark operators remain. Each type can be written as 

eabcea'b'dqa'^qaXqb'TPqb) = e a b c e a 'b>c j \5 ara 5 Vb {qT a q){qT^ q) + ^t A a ^Vb{qV a t A q){qT & q) 

+ 2 -5 a , a C(qT«q){qTH B q) + <^f b , (qT a t A q){qTH B q) j , (89) 

where T a = {7°, ^7^75} and t A is the generator of SU(3) normalized as Tr^i 5 ] = \5 AB . Combined with the product 
of cpsilon tensors e a bc£a'b'c = dbb'Saa' — dba'dab', one finds that the second and third term in the R.H.S. of Eq. 
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vanish. In the last term, the product of the generator of SU(3) can be simplified using the following identity: 

(90) 



,A ,B _ ,C , 

Vo'ii'6 _ o" T a'a l b'bT 



l a'a z b'b g° l a'a l b'b 

where only the first term in the R.H.S. of Eq. (|90|) survives after multiplying it with the epsilon tensors e a bc£a'b'c- 
Then Eq. can be simplified as follows 

e a bce a >b>c{qa'r a q a )(q b ,TPq b ) = e abc e a , b , c j ^S a , a S b , b {qT a q)(qr^q) + \t A a ,t A , {qT a t B q)(qYH B q) 

= l(qr a q)(qT f3 q) - 2(qT a t B q)(q^t B q). (91) 

One can take another Fierz rearrangement to [uT a t A u){uY a t A u) type operators in Eq. ([55)1 . Then one can obtain 
the following relations when taking the symmetric and traccless parts of the operator relations: 

(S 7 Q i A u)(u/i A u)| s , t = - ^(«7°u)(*A)l«,t - ^(u 7 a ^ 5 u)(u^ 75 u)\ s , t + l(ucr a u)(ua°Pu)\ s>t , (92) 

(u^ l5 t A u)(u^ l5 t A u)\ s , t = - ^-(uri 5 u)(u^ l5 u)\ Sit - ±(uru)(u^u)\ Sit - ±(ua a u)(ua°Pu)\ s , t , (93) 

where Sj t means symmetric and traceless. Therefore, only three independent twist-4 (dimension 6 spin-2) matrices 
remains. Using the twist-4 effects in the deep inelastic scattering data on the proton and neutron target, one can in 
principle extract two independent constraints to the three independent matrix elements. To determine all the matrix 
elements, we will additionally use one constraint adopted by Jaffe |47|: (ua a u)(ua oa u)\ s t = 0. 

2. Twist-4 operators with mixed quark flavor 

The second type of four-quark operator appearing in the nucleon sum rule are of the following mixed quark flavor 
operators 

eahcta'bc tV^tV {u T a C~i4 lv Cu al ) = e abc e a , bc , ^( 1 5 ^T k j" 1 5 )(u a ,T°u a )(d c ,T k d c ) ■ Tr [^ V CT C] 

eabcCa'bc' TS\ ~ 8q a (u al j a u a )(d cl d c ) 
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8(<7/3 7 a + 3a/3^)(«a' 7 a Wa)(4' 7 ' 9 4) 
8(9/3 7a - 3a/3^)(«a' 7 Q 75Ma)(4'7 /3 754) \, (94) 



where we have again used Fierz rearrangement to express the operators in terms of quark-anti quark type, and have 
neglected operators that are odd in parity and time reversal symmetry. 

As in the case of the pure quark flavor case, the four-quark condensates in Eq. (|94p can be decomposed into two 
different color structure according to Eq. (|9Tj) and Eq. ([90]) . We can not reduce the number of independent operators 
as in the previous subsubsection because performing a similar Fierz rearrangement as in Eq. (|92[) and Eq. (|93[) . we 
find new mixed flavor operators in (ur a d)(dr°u) type. 

3. Contributions of dimension 6 four-quark to the OPE 

In summary, the independent four-quark condensates appearing in our nucleon sum rule are given in Table |H Not 
all the matrix elements are known. 

As for the dimension 6 spin-0 (scalar) operators, we will assume the factorized form as {uu) 2 p 7 ; although this 
assumption has not been justified. Keeping only the linear density terms, they can be written as 



V K+(m q ) 



(qq)v ac ([qq\o)p P, (95) 



where '/' is a parameter introduced in Ref. \\A 
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Quark flavor 


qi = qi = q 


qi ^ qi 




(g7 Q 75g)(g7 ,3 75g)|s,t 


(qi7° , 75<7i)(g27 ,3 75<?2)|s,t 


Dimension 6 spin-2 


(qi a q)(q^q)\ s .t 


(9i7 Q <7i)(927 ,3 92)kt 


(Twist-4) 


(qa a q)(qa°Pq)\ s , t 


(9i7 Q 75i A '7i)(g27' 3 75i A 92)kt 






{qil a t A qi){q2^t A q 2 )\s,t 


Dimension 6 spin-1 




(qij a qi){q2q2) 


(Vector) 




(qi-y a t A qi)(q 2 t A q2) 




(q-y a j5q) (qj a j5q) 


(<?i7«75gi)(<?27 a 75g2) 


Dimension 6 spin-0 


{r/ a q)(qj a q) 


(gi7 a gi)(<?27 a g2) 


(Scalar) 


(qa oa q)(qa° a q) 


(gi7a75* A gi)(g27 a 75i A g2) 






(gi7 a i A gi)(g27 a t A g2) 



TABLE I: Independent four-quark operators appearing in the nucleon OPE with Iofle's interpolating current. l q%' and '52' 
represents light quark flavors. 



Dimension 6 spin-1 (vector) operator are neglected as their contribution to the nuclear symmetry energy is negligible 
to leading order in density. 

Dimension 6 spin-2 are the twist-4 operators. The twist-4 operators appearing in the nucleon sum rule have 
similar structures as those appearing in the higher twist effects in deep inelastic scattering [34| 53]- If the higher 
twist effects are measured with precision in DIS for the proton and neutron target, the nucleon expectation value of 
(u r ) a ^jc > t A u){d'y^^jr ) t A d)\ s ^ can be estimated with the same precision |33| . With further plausible arguments (Appendix 
C) on the ratio of u quark and d quark content of the proton such as those used in Eq. (|53j) . one can estimate proton 
expectation value of (u r y a t A u){u'y l3 t A u)\ s ^, (uj a j5t A u) (ujPjs^u) | S)t and (u r y a t A u)(d'y^t A d)\ s ,t- 

From these condensates, one can estimate nucleon expectation value of all the twist-4 operators for the single flavor 
case given in the first column in Table U with the extra constraint discussed above. For the mixed flavor condensates 
given in the second column, one can not deduce all the matrix elements (uj a j5u)(d r y l3 j5d)\ s ^ and (uj a u)(dj^d)\ s ^ t 
from (uj a ^f5t A u)(d^^^5t A d)\ Si t and (uj c *t A u)(d'y l3 t A d)\ Si f We will however, neglect (u-f a j5u)(d'y l3 j 5 d)\ s .t and 
(u'y a u)(d'y^d)\ s ,t in our present analysis, as these mixed quark flavor condensates do not give important contribution 
to the nuclear symmetry energy in the linear density order. 

The proton expectation value of the deducible twist-4 operators can be parameterized into the following forms: 

((9i7 Q 75^<Zi)fe7 /3 75^9 2 )) P | s ,t = (""^ - \s afi ) T U> ( 96 ) 

((gi7 a ^9i)fe7^ A <Z2)) P | s ,t = j— — -7p ( u"^ -a** ) T 2 ?2 , 



1 


M N 


4ira s 


2 


1 


M N 


4ira s 


2 


1 


Mn 


47ra s 


2 


1 


Mjv 


47ra s 


2 



u a u p - 


V') 


I T 1 


u a u^ - 




I T 2 


u a u p - 




1 T 3 

1 9192 


- 




I 9192 



<(9i7 Q 75<?i)fe7 /3 75<Z2)) P | s ,t = - \g a ^ T 3 92 , (98) 

<(gi7 a 9i)fe7%)} P kt = T^-^f (u a u? - -aA T 4 (99) 



where l q\ and '92' represent quark flavors. We have extracted the T l s from the matrix elements estimated in Ref. [331 ] 
and listed in Table HU 

Using the paramctrization of the nucleon expectation value of the twist-4 operators together with the linear density 
approximation given in Eq. (|45|) . the contributions to the correlation function from the four-quark operators can be 
written as 

2 ,_ v2 1 1 



nfL M (<Zo, \q\) = ~ (uu) 2 ^ + ___[Ti, - T 2 > 



1 


1 


M N 


<1 2 


47ra s 


2 


4 


1 


M N 


<1 2 


47ra s 


2 


4 


1 


M N 



{[T\ T\] [T\ T\]I) P -±^-^L ( [T 3 _ T%] [T 3 _ r 4 ]/} p 



ng =6 , u (^, \q\) = - ^T-^a - Tl d \p 



([To - T\\ [T\ - T\]I) p + ± J- (m - T%] [T\ T\]I) p, (100) 



g 2 47ra s 2 u u UJ L 1 11 " 3q 2 Awa s 2 
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First set (GeV 2 ) 


T 1 

' uu 


Tdd 


^UU 


Tdd 


^UU 


Tdd 


T-uu 


T 4 

J-dd 


Tud 


Tud 


Ki = Kjp 


-0.132 


-0.041 


0.154 


0.048 


0.842 


0.262 


-0.875 


-0.272 


-0.042 


0.049 


Kl = Kl d {P + l)/p 


-0.071 


-0.012 


0.070 


0.012 


0.424 


0.072 


-0.422 


-0.072 


-0.042 


0.041 




-0.042 


0.002 


0.033 


-0.002 


0.240 


-0.012 


-0.233 


0.012 


-0.042 


0.031 


Second set (GeV 2 ) 


T 1 

' UU 


Tdd 


^UU 


Tdd 


-* UU 


J-dd 


■*■ UU 


T-i4 
±dd 


Tud 


Tud 


Ku = ~Kud 


0.215 


0.124 


-0.432 


-0.265 


-1.778 


-1.091 


2.104 


1.290 


-0.042 


0.057 


Kl = -K x ud {l3 + l)/p 0.154 


0.100 


-0.337 


-0.219 


-1.336 


0.868 


1.610 


1.046 


-0.042 


0.056 


Kl = -Kit/? 


0.125 


0.085 


-0.297 


-0.202 


-1.137 


-0.773 


1.395 


0.949 


-0.042 


0.058 



TABLE II: Two sets for T's. Three different classifications of T % follow classifications of Ref. [33| . Detailed treatment is given 
in Appendix C. 

where T' = \(T l uu + T dd ) and T\ = ^(T^ u — T dd ), and all the vacuum and scalar condensates are factorized as 
Eq. ([55)1 and Eq. (p5")l . The corresponding Borel transformations are given as follows: 



Bpi,=6 ) ,(gg,|5l)] = |<uu>a fJ i« 



2 

1 Mat 
47ra s 2 



{[^ - T„ 2 d ] - ([Tq - T 2 ] - [T\ - I*]/) + i ([T 3 - T 4 ] - [T 3 ! - T^]/) | 



(101) 



sp^,^ 2 , kl)] = {pi - tU + ([Tg - r 2 ] - p* - r 2 ]/) - 1 ([r 3 - t 4 ] - F\ - t\]i)} p l~I. 

(102) 

Here, we have neglected the scaling of the matrix elements coming from the anomalous dimension of the dimension 6 
operators To„ ■ Although the twist-4 matrix elements are estimated at the separation scale of 5 GeV, and the matrix 
element we need is at lower energy scale close to the Borel mass, we will neglect the running of the matrix elements 
through the anomalous dimension for operator To n , because the present estimate of the matrix elements already 
contains ±50% uncertainty. Throughout this paper, we used a s ~ 0.5 as in Ref. |33l. |34|. 

I 

IV. NUCLEON SUM RULE AND NUCLEAR A. QCD sum rule formula 

SYMMETRY ENERGY 



The quasi-nucleon self energies in rest frame can 
be obtained in QCD sum rules by taking the ratios 
Eq. dS5])/Eq. (O and Eq. (|3D]l/Eq. for both the 

proton and neutron: 



We have expressed the self energy contributions of 
the nucleons that contribute to the nucleon energy as 
Eq. (fl~6f in terms of the Borel transformed OPE as given 
in Eqs. <J3HJ), (EH) and gO])- The next step is to substitute 
Eq. (fT6|) to Eq. ([8]) to extract the symmetry energy as de- 
fined in Eq. Then, there will be the trivial kinematic 
correction coming from the three momentum dependence 
in the kinetic energy part of Eq. (|16p . This term is uni- 
versal and corresponds to the term in Eq. ^j. Instead 
of following the full procedure, in this work, we will just 
concentrate on the contribution coming from the scalar 
and vector self energy. This corresponds to calculating 
the contribution to the nuclear symmetry energy from 
potentials in effective models. 



E, 



Ml 



J\f n 'P{ P ) 



(103) 



where subscripts q, V{I) are meant to represent the po- 
tential part of Eq. (|16[) in the asymmetric nuclear mat- 
ter. To discuss different approximations of self energies in 
terms of the density and the asymmetric factor, we intro- 
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duce the following symbols AfJ 1 '^ n\(p) and rn(p) 



m i 



MP 



2 2 



^ \P) - u {p°J°) + U ( P .I°)P + {p,I) P 



m I 



2 2 



(104) 



(105) 



where the superscript 'n', 'p' represents either the neu- 
tron or the proton respectively. For the pair of subscripts 
(p m the first index represents the order of the den- 
sity, while the the second index represents the isospin. 
Due to isospin symmetry, the iso-scalar terms have the 
following relations: 



V 



0' 



(106) 
(107) 



where I is the integer for the order of the isospin. All the 
terms appearing are summarized in the Appendix D. 

Because the dominant term of T> n ' v (p) is T>™ p p IO y one 
can expand the denominator in terms of (l/2?^o jo)) 
times condensate. After rewriting this with powers of 
p and I, one can express the potential part of a single 
nuclcon energy as 



oo oo 



fe=l i=0 



rpn,p k 
J2j V,(p k J')P 



(108) 



where E y ,v pk ^ are written in terms of N n ' p {p) and 
V n ' p (p). Averaging Eq. (fTU5|) as Eq. © and collecting 
terms of I 2 from Eq. ©, one can extract Ey m (p) as 
follows: 



;P ■ ( E v,( P ,i) ~ E v{p I) ) + -p 2 ■ (E v{p2 j } - £y )((9V) ) + jP 3 ■ ( E v,(p*,i) ~ Ey {pa j) ) + 



\P 2 ■ (Ev,(p 2 J 2 ) + • E y i(p 2 i/ 2 ) ) + -p 3 ■ (Ev^ p3 j 2) + E^ p3 l ^) + 



(109) 



For terms linear in density, one can see that the first term in the upper bracket of Eq. (|109p corresponds to the 
form given in Eq. ([9]). The explicit expression in terms of A/^ m ir , and D ? m ^ is 



E ^T = - 4 p- 



• -2A/i 



TP \ (PJ)I (J)P \ 

(p°J°) \ u (p°,I a )) 



.(- 2 V? 



(pj; 



(110) 



valid to leading order in density. 

When higher density dependence of the condensates are calculated, Eq. (|109[) provides a systematic expression of 
Ey nl (p) that includes higher p n - 2 terms. 



B. Sum rule analysis 

In principle, a physical quantity extracted from the 
QCD sum rule should not depend on the Borel param- 
eter M 2 . However, since we truncate the OPE at fi- 
nite mass dimension, such a physical quantity should be 
obtained within a reliable range of M 2 (Borcl window) 
with "plateau". While we do not find the most stable 
"plateau" with a extremum in the appropriate Borel win- 
dow, one finds that the results has only a weak depen- 
dence on M 2 . 

The well accepted Borcl window for nucleon sum rule 
is 0.8 GcV 2 < M 2 < 1.4 GcV 2 0. But as our sum 
rule contains the newly added twist-4 four-quark oper- 



ators, the Borel window needs to be re-examined. We 
determine the upper Borel window by requiring that the 
quasi-nucleon contribution to be more than 50% of the 
total sum rule so that the continuum contribution be less 
than 50%. As for the lower limit, for the same OPE, we 
restrict the contribution from the highest mass dimension 
operator to be less than 50% of the total contribution. 
For the quasi-nucleon energy in medium rest frame, we 
applied this prescription to the R.H.S of Eqs. (|38j) . ([39]) 
and (j4T)]) . 

The Borel curves for the three invariants (Eqs. ([55]) ), 
([39]) and ([30])) are plotted in Fig.Q] Here, all the graphs 
are obtained with the T l s using the K\ = K^diP + 1)//3 
estimates from the first set of Table [TTJ From the left 
figure of Fig. [TJ one can get an acceptable Borel window 
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FIG. 1: Borel window for three invariants. Left : the borel 
window for B[Tl s (qo, \q\)] and B[H u (qo, \q\)]. Right : the borel 
window for B[H q (qo, \q\)]- In both figure, the thick lines in- 
creasing with the Borel mass represent the ratio (the con- 
tribution of highest dimensional operators) /(the total OPE), 
and the thin lines decreasing with the Borel mass represent 
the ratio (the continuum contribution) /(total contribution). 
These graphs are obtained with T l s in K\ = K^diP + 
estimation from the first set of Table [II] 

for B[Tl s (ql \q])} (Eq. (J38J)) and B[Il u (ql\q\)] (Eq. gOJ). 
However, in the right figure of Fig. [TJ B{H q (q 2 , \q])] do 
not provide an acceptable Borel window. While the usual 
Borel window is obtained by requiring that the power and 
continuum correction arc both less than 50% of the total 
OPE, we will loosen the condition to be less than 75% in 
this case. 

This large power correction, may be caused by an over- 
estimated ([?9]o)vac- A s mentioned in the previous sec- 
tion, all the vacuum expectation values of four-quark 
operators are factorized as in Eq. (|95p . Only large N c 
supports factorization in the vacuum. Hence, the gener- 
alization to the nuclear medium can only be an order of 
magnitude estimate with large uncertainty. For example, 
as one can see in the right figure of Fig. [TJ the lower and 
upper boundary from B[H q (q 2 , \q\)] are already largely af- 
fected by whether the vacuum value ([<7<?]o)vac i s included 
or not. Also, the neglected twist-4 matrix elements T 3 d 
and T* d for (u^^id^d)^ and (M7 Q u)(J 7 a d)| s ,t 
may be another reason for the large uncertainty. If the 
vacuum expectation value of four-quark operators and 
T 3 d and T* d can be determined well, we can discuss 
about the stability of our sum rule in a more reliable 
way. The second set of T l s from Table QTJ do not pro- 
duce any acceptable Borel window. In conclusion, we 
will use the results from the following Borel window: 
1.0 GeV 2 < M 2 < 1.2 GeV 2 . 

In the analysis to follow, for the symmetric nuclear 
matter case, we will denote the twist-4 condensates con- 
tribution to the quasi-nucleon self energy as Et and 
the total quasi-nucleon self energy in the rest frame as 
Eq.v(i=o)- F° r the asymmetric nuclear matter case, we 
will use two sum rules for Ey m : one that includes con- 
tributions up to order p terms and another one up to 



FIG. 2: The ratios between quasi-nucleon self energies and the 
vacuum mass. The different lines represent (M^/Mn (dashed 
blue), Eu/Mjv (dot-dashed red), Y>t/Mn (dotted black) and 
Eq iV (i^o)/MN (solid black)) respectively. 

p 2 . The former sum rule will be called the linear p sum 
rule (Ey™) and the later p 2 sum rule (Ey™ 2 ). As for the 
value for the anti-nucleon pole, an optimal 'in-medium' 
value ranged -0.2 GeV < E q < -0.4 GeV will be used 
for each different estimation of twist-4 matrix elements 
in the sum rule for the quasi-nucleon self energy, while 
the 'bare' value E q = —Mm will be used in the sum rule 
for the nuclear symmetry energy. This is so because the 
quasi-hole contribution in the nuclear symmetry energy 
comes with a term proportional to the density (Eq. (JH])). 
Nuclear matter density p is set at the saturation den- 
sity po = 0.16 fm -3 and the corresponding quasi-nucleon 
three- momentum |<f| is taken to be 270 MeV, the Fermi 
momentum of normal nucleus (po = 0.16 fm~ 3 ). The 
light quark (u, d quark) mass m q is taken to be 5 MeV. 

As for the density dependence of dimension 6 spin-0 
condensates, different / values are used for every esti- 
mations of twist-4 matrix elements. For the first set 
of Table HI / = -0.2 for K x u = K\J$ (correspond- 
ing E q = -0.26 GeV), / = -0.17 (corresponding E q = 
-0.34 GeV)for K\ = K\ d ($_+ l)//3 and / = -0.12 for 
K l = K ld (corresponding E q = -0.4 GeV). This pa- 
rameter set of /s are chosen to satisfy the self consistency 
constraint as given in Eq. (|17p for the quasi-hole value. 
Again, the second set of the Table |TT] do not provide a 
set of /s which satisfies the constraint Eq. (fT7|) . Detailed 
discussion for related parameters (/ and E q ) will be given 
in later section. 



1. Symmetric nuclear matter 

First, we investigate the quasi-nucleon self energies 
in the symmetric nuclear matter with twist-4 conden- 
sates. Throughout the analysis, we check the result 
against the 1 = case. In Fig. [21 we plot the ratio 
to the nucleon mass in vacuum of the in-medium scalar 
self energy (M^/Mn), the vector self energy (E v /Mn), 
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FIG. 3: Density (right) and |g| (left) dependence of the ratios 
between quasi-nucleon self energies and the vacuum mass. 
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FIG. 4: Ey'™ (left) and its \q\ dependence (right). The unit 
of the vertical axis is GeV. 



the twist-4 condensates contribution (St/M^v) and the 
potential part of total quasi-nucleon self energy in rest 
frame (E q v(i=o)/Mis!)- For the twist-4 matrix elements, 
we take K\ = K\ d ($ + l)//3 from the first set of Ta- 
ble |IT] and corresponding E q = —0.34 GeV, which gives 
the average result. From our analysis, we find that the 
contribution of the the twist-4 condensates give enhance- 
ment of the quasi-nucleon self energy by ~ 50 McV. 
When / = 0, we find the ratio E q ^v(i=o)/^N — 0.96, 
M* N /M N ~ 0.62 and Y> V /M N ~ 0.33. By using afore- 
mentioned parameter set for / < and E q , the ratios 
become E qy{I=0) /M N ~ 0.87, M^/M N ~ 0.56 and 
Y, v /Mn ~ 0.30, which are comparable with previous 
studies [TMl. When the second set of Table HD is used 
for the twist-4 matrix elements,we do not find a stable 
behavior in the same Borel window in contrast to the 
case with the first set as shown previously. By setting 
/ > for the second set, E q y^ I= ^/M^ can be adjusted 
to ~ 0.9 which is a typically acceptable value. But even 
so, there is no reasonable / and E q for M^/Mn and 
S„/Mjv which satisfies Eq. (jTTJ) . The estimates for T's 
given in the second set of Table [XT] do not reproduce the 
aspect of the nucleon sum rule that is consistent with 
the Driac Phenomenology [ll|. Hence, we will continue 



the present analysis with estimates for T l s given by the 
the first set of Table HJ 

The quasi-nucleon three momentum dependence is 
plotted in the left figure of Fig. [3] for / < cases: 
one finds that the ratio E^/Mjv and St/M/v do not de- 
pend strongly on the quasi-nucleon three momentum. On 
the other hand, M^/Mn shows significant change when 
\q\ > 0.5 GeV as in Ref. So this sum rule analysis 

works in < \q\ < 0.5 GeV region, which is consistent 
with our phenomcnological ansatz that assumes a mo- 
mentum independent self energy. 

As all the condensates in our nucleon sum rule are 
estimated up to linear order in density, the results may be 
valid at least near the nuclear saturation density region. 
In the right figure of Fig. [3J the density dependence of 
the quasi-nucleon self energies is plotted for 0.4 < p/po < 



1.6. Here we used the parameter set / = —0.17 and 
E q = —0.34 GeV determined at the saturation density, 
as our nucleon sum rule do not depend strongly on E q 
in —0.6 GeV < E q < —0.3 GeV which covers naively 
estimated E q in 0.4 < p/p < 1.6. Yj v /Mn and T,t/Mn 
are enhanced as the nuclear matter becomes denser but 
M* N /M N is reduced. 



2. Asymmetric nuclear matter 

In our sum rule, nuclear bulk properties in the asym- 
metric nuclear matter are parameterized by the asym- 
metry factor /. If one plots the quasi-nucleon self energy 
as a function of / to leading order in density, Ey™ can 
be obtained from the difference between the slope of the 
quasi-neutron and the quasi-proton (Eq. (|110p ). 

Ey™ is plotted in Fig. |H One finds that Ey ™ ranges 
from 40 to 50 MeV, which agrees with previous studies. 
The results in the figure also show that including the 
twist-4 contribution slightly enhances the nuclear sym- 
metry energy. 

In the right figure of Fig. [U one finds that Ey™ do not 
depend strongly on the quasi-nucleon three momentum 
up to 0.5 GeV. This result agrees with the quasi-nucleon 
three momentum dependence of the quasi-nucleon self 
energy. When the second set of Table [H] are used for the 
T's, we find that -Ey ™ depends strongly on the quasi- 
nucleon three momentum than the first set. 

One can also work out Ey 2 , although with larger un- 
certainty than that of Ey™. The density dependence of 
£^™ 2 and E% m in 0.4 < p/p < 1.6 are plotted in Fig. [5] 
Here again, the four-quark condensates give a important 
change to the density behavior of Ey m . For / = 0, the 
contribution of T l s gives slight enhancement of Ey ™ 2 at 
higher nuclear matter density but for / = —0.17, it gives 
reduction of Ey ™ 2 at higher density. This means, in the 
symmetry energy case also, the scalar four-quark opera- 
tors in our sum rule may play a important role in reduc- 
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FIG. 5: Density dependence of E aym . Left : density depen- 
dence of Ey™2. Right : density dependence of _Ej^ m . The 
unit of the vertical axis is GeV. 



FIG. 7: Sensitivity analysis with ([gsq^u ■ Gq]o)p- Left : 
{[g a q*o ■ Gq]o) P dependence of E q , v{1 = 0) . Right : ([g a q*o • 
Qq]o) P dependence of Ey ™. The unit of the vertical axis for 
the right figure is GeV. 
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FIG. 6: Scalar- vector self energy decomposition of Ey™. Left 
: Ey™ without twist-4 contribution. Right : Ey m with twist- 
4 contribution. The unit of the vertical axis is GeV. 



tion mechanism of the nuclear symmetry energy. How- 
ever, E^ m is slightly reduced by T l s as twist-4 matrix el- 
ements enhances M^/Mm- The parameter set / < con- 
tributes differently Ey™ 1 and E s ^ m ; reduction for Ey™ 1 
and enhancement for E s ^ m in 0.4 < p/po < 1.6. 

In Fig- El wc plot the scalar [n^'j^j] ) anc ^ vec tor 



In the left fig- 



\B[U q (ql,\q\)\) Sel1 eiier gy P ait OI ■ 

ure, we plot the result without the twist-4 contribution, 
while in the right hand figure, we include the contribution 
from twist-4 matrix elements. While both the scalar and 
vector self energy give positive contribution to the self 
energy in the left hand figure, one finds that in the right 
hand side, the scalar and vector give negative and posi- 
tive contributions respectively. The right figure is consis- 
tent with the general trends in RMFT results [2!| , which 
shows scalar self energy part gives negative contribution 
and vector self energy part gives positive contribution 
from the exchange of 8 and p meson exchanges respec- 
tively. One can infer from this result that the twist-4 
contribution mimics the exchange of the 5 and p meson 



and that it constitutes an essential part in the origin of 
the nuclear symmetry energy from QOD. 



3. Uncertainties 

In general, there are two quark-gluon mixed opera- 
tors with contracted spin indices, ([g s q<J ■ Gq]o) P and 
([g s q*<j ■ Gq]o) P , which are not accurately determined. 
([g s q&'Gq\o)p does not appear in our sum rule with Ioffc's 
nucleon interpolating field (Eq. (JTHJ) ) . As for the oper- 
ator (\g s q^a ■ Gq]o) P , the proton expectation value has 
been estimated in Refs. [lj, [3 5H, [4(| to be in the range 
of -0.33 GeV 2 < ([^gV • Gq\o) P < 0.66 GeV 2 . Hence, 
we investigate the ([g s q^cr ■ Gq]o)p dependence in Fig. [7] 
The matrix element ([g s q^a ■ Gq]o)p do not give impor- 
tant contribution to the quasi-nucleon self energies in the 
range -0.33 GeV 2 < ([g s q^a ■ Gq]o) P < 0.66 GeV 2 as we 
arc not interested in the accuracy of 10 McV. However, 
for Ey™, ([gsq^a ■ Gq]o)p such magnitude gives nontriv- 
ial fractional change to Ey 1 ™. We choose the value as 



{[g s q^o- -Gg 
in Ref. [H, 



-0.33 GeV 2 in this study as was done 



Conclusion 



In this paper we studied the nuclear bulk properties in 
the asymmetric nuclear matter by calculating the quasi- 
nucleon self energies with QCD sum rule approach. In 
particular, we identified all the twist-4 local condensates 
appearing in the nucleon sum rule. Using the existing es- 
timates for the twist-4 matrix element from DIS, we were 
able to find the magnitudes of all the twist-4 matrix ele- 
ments (T l ) in our sum rule except for two mixed quark 
flavor type condensates. We have calculated the nuclear 
symmetry energy and found that twist-4 contributions 
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are non-negligible and essential to give a phenomcnolog- 
ically consistent result with RMFT for the quasi-nucleon 
self energy and the nuclear symmetry energy. 

For the symmetric nuclear matter case, we found that 
E q y[i = Q) is enhanced by ~ 50 MeV with T z s in the first 
set of Table [TTJ As the T l s in the first set of Table M 
provides qualitatively reliable sum rule results while T l s 
in the second set of Table [TT| do not, we conclude that 
taking the sum rule results with T l s in the first set is 
reasonable choice. With parameter set / < 0, dimension 
6 spin-0 (scalar) operators reduces E qi vn = ff\/M^r to ~ 
0.87. 

For the asymmetric nuclear matter case, we confirmed 
two meaningful facts. First, QCD sum rule technique can 
be used to successfully reproduce the acceptable result 
for the nuclear symmetry energy at the nuclear matter 
density. Second, dimension 6 spin-2 (twist-4) conden- 
sates plays an essential role in making the scalar part 
contribute negatively to the self energy and thus provid- 
ing a consistent picture for the E sym with the RMFT 
results HI: 



risym _ 1 



f P - fs 



El 



PB, 



(111) 



where f p is isovector p meson coupling, fg is isoscalar 
5 (/o) coupling and is the nuclear matter density. 
Moreover, our approach provides a first attempt to un- 
derstanding the origin of Ey™ in terms of local operators 
directly from QCD. This extends the analogy between 
QCD sum rules to RMFT for the symmetric nuclear mat- 
ter established in Refs. [ill EH Ej], to the asymmetric 
limit. 

While the uncertainties in T z s and four-quark scalar 
with parametrization / are still large, attempts to mea- 
sure the twist-4 contribution in DIS at the future upgrade 
at Jefferson Lab is expected to lower the uncertainties 
and provide more insights to the value for the nuclear 
expectation value of four-quark operators. 
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Appendix A: Baryon octet mass relation 

In this section, we summarize an essential argument for obtaining ([qq]i) p = \ ({p\uu\p) — {p\dd\p)^j from Ref. [22[ 
and Ref. [1^] . Phenomenologically, the nucleon mass can be expressed in terms of the matrix element of the trace of 
the energy-momentum tensor: 

mjf(N\M N \N) ={N\6^\N). (Al) 
Using the equations of motion, the trace of the energy-momentum tensor can be written as 

O^p = m u uu + m d dd + m s ss + myjih + ■•• 

- — ) G 2 + m u uu + m d dd + m s ss + 0(p 2 /4m?), (A2) 
4a s / 

where h are heavy quark fields, and the gluonic term comes from the trace anomaly [4lTj43l |. j3 = —9a^./2ir is the 
"reduced" Gellmann-Low function in which heavy quark contribution has been subtracted using the heavy quark 
expansion |44j . 

Eq. (|A2j) can be applied to the lower-lying baryon octet. The lower-lying baryon octet mass relation with SU(3) 
flavor symmetry is as follows: 

m p = A + m u B u + m d B d + m s B s , 

m„ = A + m u B d + m d B u + m s B s , 
m s + = A + m u B u + m d B s + m s B d , 
m s - = A + m u B s + m d B u + m s B d , 
m H o = A + m u B d + m d B s + m s B u , 

to s - = A + m u B s + m d B d + m s B u , (A3) 

where A = ((B /4a,,)G • )„, B u = (uu) p , B d = (dd) p and B s = (ss) p . In this relation, correction terms for hyperon is 
neglected [28| . From (|A3|) one can obtain 

(p uup) - {p ddp) = , \^ , ■ A4 

2m s - (m„ + m d ) 
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Appendix B: A simple constraint for twist-4 operators from zero identity 

In this section, we show an explicit calculation for a simple constraint using the zero identity (35j . For pure quark 
flavor diquark structure, 

e abc {u T a CTu b ) =0 If (CT) T = -CT, (Bl) 

where T is Clifford basis. So (T = {/, 75, 27^75}) satisfies the above condition. Therefore, a constraint for four-quark 
operator can be obtained by requiring that the Fierz transformed form of the products of above diquarks are zero. 
An example is the following: 

eabc£a>b>c(UaCj^5U b )(uv"/ u "/ 5 Cu a >) = tabc^a'b'c — {u a >T° U a )(u b 'T k U b ) ■ Tr ["/^T k J v CT^ C] 

= e abc e a ' b ' c ^(u a 'U a )(ub'U b ) ■ (4^) 

+ (u a >j 5 u a )(u b ,j 5 u b ) ■ (-4g M1/ ) 
+ (u a >j a u a )(u b ij p u b ) ■ (AS^p^a) 

- {U a >i a i5 u a){u b >"f r: 'j5Ub) ■ (^S^pua) 

- (u a ><j aa u a )(u b ,<T f:S Pu b ) ■ -Tr [l^crp^ u <j a& ] 

+ (u a >i a U a )(u b >-f f:S J 5 U b ) • (8^ 

+ (u a >u a )(u b -(T aa u b ) ■ (8ig afl g a u) 

+ {u a >l5Ua)(u b >cr a& u b ) ■ (4e^„ aai ) j = 0, (B2) 

where S/j, al/ p = g^ a 9vp + g^f)9av — 9^9a/3- By subtracting Eq. (|B2[) from Eq. (|87[). Eq. (|57|) can be simplified into 
Eq. (USD- 



Appendix C: Estimation of twist-4 matrix elements 

In this section, we provide a detailed treatment for extracting T l s from the values estimated in Ref. [33| . In Ref. [33j |. 
twist-4 operators which appear in our nucleon sum rule are given as 

1 M N ( 1 



47ra s 2 V 4 
1 M N ( 1 



a up - ^g a [j ) Kl = ((u^ a -f 5 t A u)(uj^ 5 t A u))p\ at t + ((u7a75* A- «)(d7/375^d)}p|s,t, (CI) 



ATra s 2 
1 M N I 



[uaUfi - ^gafs^j K\ = ({w-f a t A u)(ujpt A u)) p \ Syt + ((wy a t A u) (d^ pt A d)) p | S)t , (C2) 



r (u a up - ^9 a pj Kl d = 2{{u lalb t A u){d"if^ b t A d)) p \ Stt , (C3) 



47ra s 2 

where we changed the normalization for the nucleon state appearing in Ref. [33| to the following: 

(N(p)\N(p')) = ^-(2*)H*{p-pi), (C4) 



with Up = po = yjp 2 + M'^. Here only K* d is uniquely determined: K^ d = —0.083 GeV 2 . One can set a constraint 
\K d \ = \K^\ ■ (3 < \K^ d \ < \K*\ with an ansatz that the ratio K l d /K l u is equal to the momentum fraction of the d and 
u quarks in the nucleon: 

J x(d(x) + d(x))dx _ 
Kd/K - J x(u(x) + u{x))dx - P ~ °' 476 - (C5) 

Varying K\ in the constraint above, one can estimate K\ as a functions of K\ (Table [TTT|) . 
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K 2 


Ki 




Kl 


KJP = -0.173 


0.203 


K 1 - 


0.083 


-0.181 


Kl d {P + l)/2p = 


-0.112 0.110 


-KdW 


+ l)/2/9 = 0.112 


-0.225 


K\ d = -0.083 


0.066 




= 0.173 


-0.318 



(GeV 2 



TABLE III: Tabic for Kl from Ref. W 



T* u and T dd can be easily estimated by taking T^ d = \K\ d from K\ and K d : 



T 1 



Tdd — Kd ^udi 

TL = \Kd = -0-042 GeV 2 



(C6) 
(C7) 

(C8) 



Similarly, one can try to obtain T 2 U and T dd from if 2 and K d . As T 2 d = \K^ d has not been determined uniquely 
as ^ud ~ \^udi we assume d that the ratio T^ u /T dd is equal to T^ u /T dd . Then, by the following relation, one can 
estimate T qq s: 



rp2 

1 dd 



{Kl - K\) ■ 1 



rpl 
± dd 

rpl 

1 uu 



rp\ ] ' UU 1 

UU / 



1 



Tld = 2 ( [A '« + K * ] ~ + T *ti 



(C9) 
(C10) 
(Gil) 



where if 2 — K d and K 2 + K d can be obtained from Table IIHI 

For pure quark flavor case, T qq and T qq can be obtained from Eq. (j9"2"]l and Eq. As discussed in Ref. [I?}, we 

neglect {ua a u){ua oa u)\^ t . Then, and T g 4 g can be related as 



rpS 
11 



rp4 

QQ 



15 i 9 2 

4 J <?<? 4^99 



15, 



-T 



</</ 



V. 

4 w 



(C12) 
(C13) 



X^s can be classified as three different classification of K l s in Tabic ILTfl K\ = {K\ d j '(3, K\ d {fi + l)/2/3, K^ d } and 
if„ = {— K* d , — K^ d ((3 + 1) /2/3, — K^ d /f3}. T qq s are classified in Table ILT1 according to these three classifications in 
two sets. 



Appendix D: QCD sum rule formulas for an< 3 Ey'™' 

In this section, we provide the detailed description for 



E, 



1,V{I) 









1 


U (P°,I°) 






I 



2Af[ pI) ) - 



TfP 1 (P' 1 ) J (TP \ 

_(p°J°) \ u (p a ,I°)> 

in QCD sum rule formula. In this formula, A/"^ and are as follows: 



(pj; 



(Dl) 
(D2) 



21 



^ p oja } = - ■^ 5 (M 2 ) 2 E 1 ([qq] ) yaa , (D3) 

A^ >J0) = - ■^{M 2 ) 2 E l {[qq] ) p - ^<?{[q{iD a iD }q] ) p L-* 

+ 3^2 {M^dqU^pErL-l + ^<f ([q{ 7o iD iD }q] ) p L-i - -^M 2 ([g s Ja ■ g q ] ) P E L-* 
+ MlJ- ^M 2 ([q{^D } q ] a ) p E L-i - JLa/ 2 ^ [{u . of + {u . G) 2 ]) p E L-* 

[TL - T 2 d ] + [Tq - T\] - i[T 3 - T 4 ]) L-*|, 



1 Mjv 



A/? pJ) = - ^{M 2 ) 2 E 1 {[qq] 1 ) p - A^ ([g - {liW}(Z ] l)pL -f 

r ([g t g ] 1 ) p S 1 L-* - ^f([q{ 7o iD iD }q\ 1 ) p L-i + -l-M 2 {\g s Ja ■ Qq]i} p E L~ 



(D4) 



P (p,/o) = - (J^2 M2E ~ ^f) ([q{^Do}q}o) p L~i 

+ *h M2 < v G2 ) P + ^ - O • G > 2 + (u • G ^> P L ^ 
+ f / • (gg)vac(fe]o> P i* - ^ ^ (p^ - ^ - m - t%] + \m ni) 

+ E q ^M 2 E L-*([ q U]o) P - (l - ([g{ 7 o^oi£»o}?]o>p,/i-* 

2 1 4 1 

- ^2 ([g{7o«A)i.Do}gM P ^ + i8^<[ff«? t<7 ' Glh)pL-^ |, (D7) 

P (p,/) = ^2M 2 Eo([q{^oiD }q] 1 ) p , I L-i 



4 7Z-(m„) , ,. , , 4 1 Af)v /r 1 9, l r 

3' ■ ^<^»>^ - 1^ - T2 1 - 3^ - ) 

£ 9 {- ( [g{ 7 oiA)iA> }g]i>p£ _ * + 3^2 <[»«« tCT ' $ u V)p L ~ f ' 
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